rable temperatures, the crystallization of babmgtomte requires 
more hydrous conditions , lower C0 2 , and slightly higher 0 2 gaci 
ties in the fluid phase than ilvaite. Since similar temperatures, C0 2 
pressures, and oxygen fugacities induced within skam deposits exist 
on Venus, ilvaite and perhaps babingtonile could have also formed 
on the surface of this planet by the interaction of the venusian 
atmosphere with extruded basaltic rocks. One factor that mig t 
mitigate against the formation of these calcic Fe 3+ -Fe silicates on 
Venus, however, are the high abundances of Mg and A1 measured 
during the Venera 13/14 [32] and Vega 2 [33] missions. The Mg 
and Al 3 * cations are not accepted into the crystal structures of ilvaite 

and babingtonite. , 

Discussion: Although magnetite is generally regarded to be 

the predominant ferric-bearing mineral on Venus, other muted- 
valence Fe^-Fe 3 * minerals known to exist on the surface of Earth 
could be stable in the venusian atmosphere. Thus, in addition to 

laihunite(whichis probably mctastable) and ilvaite and babmgtomte 

(both of which may be found in rocks depleted of Mg and Al), oxy- 
amphiboles and oxy-micas may also be major constituents of the 

venusian surface. The opacities and high electrical conductivities of 

such mixed-valence Fe^-Fe 3 * silicate minerals, the properties of 
which resemble magnetite [34], may also contribute to high radar- 
reflectivity regions in the highlands of Venus [35]. 
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THERMAL BUOYANCY ON VENUS: PRELIMINARY RE- 
SULTS OF FINITE ELEMENT MODELING. J. D. Burt and 
J. W. Head, Department of Geological Sciences, Brown University, 
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Introduction: Enhanced surface temperatures and a thinner 
lithosphere on Venus relative to Earth have been cited as leading to 
increased lithospheric buoyancy. This would limit [ 1 ] or prevent [2] 
subduction on Venus and favor the construction of thickened crust 
through underthrusting. Underthrusting may contribute to the for- 
mation of a number of features on Venus. For example, Freyja 

Montes, alinearmoun tain belt in the northern hemisphere of Venus, 

has been interpreted to be an orogenic belt [3] and a zone of 
convergence and underthrusting of the north polar plains beneam 
Ishtar Terra, with consequent crustal thickening [4]. Such mountain 
belts He adjacent to regions of tessera and contain evidence of 
volcanic activity. Tessera is also considered to consist of thickened 
crust [5] and crustal underthrusting is one possible mode for its 
formation [4]. Models for the formation of the mountain belts and 
associated features must then explain compressional deformation 
and crustal thickening, as well as melt production. 

In order to evaluate the conditions distinguishing between 
underthrusting and subduction, we have modeled the thermal and 
buoyancy consequences of the subduction end member. This study 
considers the fate of a slab from the lime it starts to subduct, but 
bypasses the question of subduction initiation. Thermal changes in 
slabs subducting into a mantle having a range of initial geotherms 
are used to predict density changes and thus their overall buoyancy. 
Finite element modeling is then applied in a first approximation of 
the assessment of the relative rales of subduction as compared to the 
buoyant rise of the slab through a viscous mantle. 

Subduction Model: In the mode l, slabs, having a thickness set 

by 90% of the basalt solidus, subduct at a 45° angle into the mantle. 
The initial geotherms match surface thermal gradients of l0°C/km, 
l5°C/km, and 25°C/km [6]. Slabs heat via conduction, crustal 
radioactivity, phase changes, and adiabatic compression. Phase 
changes involving the conversion of basalt to eclogite at depths of 
60 to 160 km and then enstatite to forsterite plus siishovite between 
260 and 360 km generate 0. 1 3 x 1 O' 5 crgs/cm 3 s and 0.36 x 1 (h 5 ergs / 
cm 3 s respectively [7] . The slab radiogenic heal production is 2 .63 x 
10~ 7 ergs/g s [8]. Adiabatic compression adds 0.5°C per kilometer 
of depth. Convergence rates ranged from 5 mm/yr to 100 mm/yr. 

The thermal evolution of the slab is followed using a finite 
difference technique [7,9]. The model region measures 800 km 
horizontally by 400 km deep. Processing ends when the slab tips 
reach a 300-km depth, implying time intervals ofl 0 m.y . to 1 00 m.y. 
Slab density changes derive from the thermal results through 
calculation of the thermal expansion (a^, = 3 x lCh 5 / 0 K [8]) and the 
effects of pressure (b = 1 x KH/kbar [8]) on initial densities set for 
zero pressure and temperature. The assumed initial density structure 
includes a 10-km or 25-km basaltic crust (density = 3.0 g/cm 3 ), a 
corresponding 25-km or 65-km-thick depleted mantle zone (den- 
sity = 3.295 g/cm 3 ). and an underlying undepleted mantle (den- 
sity = 3 .36 g/cm 3 ). Density changes due to the phase transitions are 
also included. Results take the form of density distributions within 
the model region. 

Finite element modeling is then employed to gauge the rate at 
which a slab having the density structure derived above will move 
through a viscous mantle. Buoyancy body forces are applied to a 
slab having a viscosity of 10 11 Pa s surrounded by a mantle with a 
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viscosity of 10'’ Pa s. Zwo stress boundary conditions are applied 
A sides of the model region, while no motion perpendrcular to 
rrpion boundaries are allowed except at the top. 

Resets: Figure 1 shows a typical result of the computations 

modeling the density changes. Densitycontours (0d ^ 

clearly delineate the slab and its crustal layer. In this case (10 Ukm 

“ fli, L™di»g n~fl.Bdo.tb. 

slabs will rise through the mantle at a rate of 5 to 10 Vm/m.y.Th 
analysis considers only the instantaneous velocity of the s'ab an 
£ not incorporate the full results of die density modeling or the 

bv negatively buoyant slabs or hindered by slabs that are positively 
buoy L Positive net buoyancy is found above die basalt-ec ogi e 
nhJe change, tending to oppose subduction. Negative net slab 

buoyancy for the hiU-lengthslab was found for allcond.l.onsw^ 

neutral buoyancy was achieved for slabs at alengthof about275 km. 
S the si ab must penetrate deeply intothemandebeforenegat.ve 

« help £*e subduction. The rate of 0. slab's ta»~ 
rise Lough the mantle isthen important in determining whether die 
slab may descend deeply enough to become negatively buoyan 
Preliminary results of finite element modeling indicate the slab 
mavrise at rates between 5 and lOkm/m.y. Thus, subducung slabs 
will tend to rise into an underthrusting position ,f *etr ’ubducuon 
rate is slow However, it may be that moderate to high rates of 
subduction will overwhelm the buoyant rise of a slab. This could 
lead to slabs being forced through the basaU-eclog.tc phase irans^ 
lion and to great enough depths to become negatively buoyant, th 
oossibly producing a self-sustaining subduction system. 

Thele initial results must be considered in light of the presump- 
tionof subduction made in undertaking the analysis. Some process 
still must be found that would carry the slab downward desp.ie.ls 
initial positive buoyancy. Further work will model more closely Je 
dynamics of the subductive motion of the slab and the e ec o 
slab density evolution on slab buoyancy, its rate of nse through the 

mantle, and the continuance of subduction. 

These results indicate that for all cases of assumed Venus 
e eo therm a lithospheric slab whose subduction h “ bccn ,nltl " = 
till instead be forced to underthrust the overriding lithosphere if the 


subduction rate is slow. This could then lead to crustd du^enmg^ 
melting, and volcanism. and possibly prov.de one models exp 
the association of compressional mountain belts ^blodcsofh.gh^ 
standing tessera, with apparent flexural rises and foredeeps, and 
with large volumes of volcanic deposits. 
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Fig. 1. Final density distribulion for 10°C/km geolherm. 25 km crust, and 5 
km/m y. subduction rale. Contours have a 0.1 g/cm spacing. 


Lava channels are a common feature in the volcanic regions of 
ihe Moon, and have now been observed on Venus ( 1 ] . TTtere has be n 
much debate about the origin of lunar channels: Are ey e res 
of erosional (either thermal or mechanical) or constructional p 
Z ^ It is necessary to determine the criteria to d.sunguish 
between the different types of channels. The clearest evidence ,s 
2 presence of levees indicates that the channel experienced a 

Rifle. » A* M». w- 

is its location along the crest of a ridge. In addition. Hadley R.llemd 
other lunar mare sinuous rilles are discontinuous, suggesung lhat 
their origin was. in part, a lava tube that has subsequently undergone 
partial roof collapse. Carr [3) and Head and Wilson [41 have argued 
Lt these rilles were produced by lava erosion. For lunar highland 
channels, which tend to be larger than their mare counterparts, 
mechanical erosion of the megaregolith is a possible process. 

Channels of several different types have been observed on the 
surface of Venus [1 ). They are probably formed by more than one 
process. They range in size from a fewkilomelersto over 00km 
[ 1 ]. The relatively short C‘tadpolel.ke") channels [5] (e.g.. 24 S 34 /) 
appear similarto lunar mare sinuous rilles in mo^^ogy. ^ey ^e 
so like certain constructional teircstrial channels (e.g. Kalaup.pa, 
Hawaii [6]) that it appears reasonable to say that they too are 
constructional channels or collapsed lava tube systems^ 

However, the longsinuous channels referred toby Baker et ah [ 1 
as “canali" pose a different problem in the understanding of then 
formation. One example of a channel of this type in the southeast 
region of Aphrodite Terra appears to show both er ° s ‘°"£ T* 
constructional characteristics. This channel is re P rM « ,,c * " ' 8 ’ j 

It is approximately 700 km long with an average width of about t 1 
km It drops a distance of 700 m from beginning to end/which means 
that die average slope is 0.06”. Its source may have been a graben 
situated at the northwest end of the channel. It appears to have 

different origins along its length. , , 

The lack of levees near the source suggests that the channel 
is erosional in this region. An inferred profile is shown as AA m 

Fig. 1. 


